Thermodynamic formalism for dissipative quantum walks 
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We consider the dynamical properties of dissipative continuous time quantum walks on directed graphs. 
Using a large deviation approach we construct a thermodynamic formalism allowing us to define a dynamical 
order parameter, and to identify transitions between dynamical regimes. For a particular class of dissipative 
quantum walks we propose a new quantum generalization of the the classical pagerank vector, used to rank the 
importance of nodes in a directed graph. We also provide an example where one can characterize the dynamical 
transition from an effective classical random walk to a dissipative quantum walk as a thermodynamic crossover 
between distinct dynamical regimes. 
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I. INTRODUCTION 

For a system in thermodynamic equilibrium Statistical Me- 
chanics and Thermodynamics provide the formalism to char- 
acterize different phases of matter, and the transitions from 
one phase to another [ 1 ] . Quantum fluctuations alone can give 
origin, at zero temperature, to quantum phase transitions [2] 
whose effects can experimentally be observed in a critical re- 
gion at small but finite temperatures. Away from thermody- 
namic equilibrium engineered open quantum systems provide 
a more recent source of interest tfH-H]]. In these latter settings 
a controlled environment can be used to drive the system to a 
steady state whose properties depend on system and environ- 
ment parameters. This approach has been used to study dissi- 
pative phase transitions controlled by the coupling constants 
in the equation of motion and to formulate a dissipative 
model of quantum computation JH . 

Most of the research on quantum dissipative systems fo- 
cuses on the static properties of the steady state, in this work 
we adopt a different dynamical perspective, where the time 
parameter plays a more explicit role in characterizing the sys- 
tem. This allows us to study the occurrence of dynamical 
crossovers or dynamical phase transitions driven by fluctua- 
tions in the ensemble of possible evolutions. The model sys- 
tems that we consider, of relevance both to quantum compu- 
tation and condensed matter, are dissipative Continuous Time 
Quantum Walks (CTQWs). The standard CTQWs (the uni- 
tary analog of classical random walks) describe the coherent 
motion of an excitation on the vertices of an undirected graph 
fffll- They are a powerful tool in quantum computation, 
providing a useful approach to algorithms' design, and an al- 
ternative framework for universal quantum computation ifTTIl . 
More recently CTQWs have also been used to describe trans- 
port properties in complex networks lll2l IT3I . while dissipa- 
tive CTQWs have found applications in the context of quan- 
tum biologv lfl4l-[l9|l. a nd quantum state transfer in supercon- 
ducting qubits H20L 12111 . 

In this work, following [22], we will refer to dissipative 
CTQWs as Quantum Stochastic Walks (QSWs). Using the 
quantum jump approach [23] we study the statistical proper- 
ties of the ensemble of trajectories originated by QSWs. The 
mathematical tool allowing us to probe the dynamics of the 
system and the structure of the space of trajectories is known 
as Large Deviation (LD) theory UM . The latter finds applica- 
tions in different areas of statistical mechanics and computer 



science (see B24H and references therein). Using the LD ap- 
proach it is possible to construct a formalism, for the ensemble 
of possible evolutions, resembling standard thermodynamics. 
Since the LD theory overlaps with the thermodynamic for- 
malism developed by Ruelle and others 1I25I1 in the context 
of chaos theory and dynamical systems, in the following we 
will use the expressions thermodynamic formalism and LD 
approach interchangeably. 

From the application of the thermodynamic formalism to 
QSWs we obtain two main results: first, a quantum general- 
ization of the pagerank vector [26, 27], a measure used to rank 
the importance of nodes in a graph; second, a tool to study 
transitions from classical to quantum dissipative dynamical 
regimes from a thermodynamic point of view (see 112 811 for 
a different approach). In particular, we provide an instance of 
a simple QS W presenting a dynamical crossover between the 
two regimes. 

The rest of the paper is divided as follows: in the second 
section we define the notion of QSW as introduced in B22I1 : 
in the third section we provide a brief introduction to the LD 
approach; in the fourth section we apply the thermodynamic 
formalism to describe the dynamics of a QSW on a directed 
graphs; while the last section is devoted to discussions and 
conclusions. 



II. QUANTUM STOCHASTIC WALKS 

QSWs are a class of stochastic processes introduced in [22] 
with the purpose of generalizing both continuous time quan- 
tum walks, and classical random walks 112911 into a single 
framework. A QSW is specified by a quantum master equa- 
tion in Lindblad form [30] describing the evolution of the sys- 
tem's density matrix 



Nl 1 

i[H,p\+Y,L kP Ll--{L{L k ,p}, 



(1) 



k=l 



where [•, •] and {•, •} stand respectively for the commutator 
and the anticommutator, and the operators L/. are called jump 
operators. The dynamics described by eq. is character- 
ized by the interplay between a coherent part associated to H , 
and an incoherent part associated to the jump operators {-£/&}. 
If only the coherent term was present the QSW would cor- 
respond to a standard CTQW. On the other hand, it can be 
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shown that if only the incoherent term was present the QSW 
would correspond to a classical random walk for the diagonal 
elements of p(t) 1E2I1 . 

A standard stochastic method to study the dynamics de- 
scribed by eq. ^ is provided by the quantum jump approach 
J23L [30ll . In this framework one considers an ensemble of 
piecewise deterministic evolutions of pure states {\ip a (t))} 
(where a labels an element in the ensemble). The average 
over different trajectories (or histories) provides the state of 
the system at each instant of time 

pft) = [i^ft))(^ft)i] ave ^E ^T a(t) ' (2) 

a— 1 

The dynamics of each pure state evolves deterministically for 
some time-interval according to the following non-hermitian 
effective Hamiltonian 
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H eB = H-iJ2~Ll L k 



(3) 



fe=i 



The deterministic evolution is interrupted at certain random 
times t by a quantum jump operator Lk, which projects the 
state according to 



\il>(r)) 



\\L k mr))\\ 



(4) 



In between random jumps the state evolves according to H e s. 
The use of an ensemble of quantum trajectories associated to 
the QSW in eq. (0) is also referred to as the stochastic unrav- 
eling of the quantum master equation. More details on this 
method can be found in I23I [30ll . 

Since the quantum walker evolves on a graph we also need 
to specify how the graph's structure enters into eq. (0). A nat- 
ural way to proceed is to use the set of multi-indices 
(where i and j are the nodes of the graph) to label the index k 
in L k 



L k -> Lij = Mij\i)(j\, 



(5) 



and the matrix M contains the rate coefficients of the jumps 
from \j) to \i), in the site basis. This way the quantum jump 
in eq. (|j) corresponds to hopping from one site to another at 
a rate given by M. For the coherent part of the evolution 
the Hamiltonian H in eq. (Q]) is identified with the adjacency 
matrix of the graph with all directions removed, so that H is 
a symmetric matrix. Note that only the dissipative part in the 
master equation is the one that encodes for the direction of the 
edges. In the next section we introduce the thermodynamic 
formalism used to study the evolution of a QSW. 



III. LARGE DEVIATION THEORY AND 
THERMODYNAMIC FORMALISM 

The dynamics of a stochastic system provides many use- 
ful information that cannot be accessed from the knowledge 
of the stationary state alone. The large deviation approach is 



the mathematical framework allowing us to access this extra 
information [24j]. In fact, recently, its application to classical 
and quantum nonequilibrium dynamics has being crucial to 
show the presence of dynamical phase transitions and dynam- 
ical crossover phenomena Il3ll433ll . In particular, as first ob- 
served in 1I32I1 . the quantum jump approach provides a quan- 
tum stochastic process that is well suited to be analysed with 
the LD theory. In what follows we show how to adapt the 
method developed in [32] to the context of QSWs. 

Given a quantum stochastic walker evolving according to 
eq. ([l]), we are interested in the time record of jump events 
associated to one or more operators L k - In the quantumjump 
approach this record is one particular realization of a quantum 
trajectory. Introducing a variable K, which counts the number 
of events (or jumps) up to a certain time for a particular jump 
operator, we define the probability 



p(K,t)=Ti[p(K,t)}, 



(6) 



corresponding to the observation of a number K of events up 
to time t. p(K, t) is the density matrix of the system condi- 
tioned on the observation of those K events. The probability 
generating function of the distribution p(K, t) is given by 

00 00 
Z(s, t)=J2 e~ sK P(K, t) = TrJ2 e- sK p(K, t). (7) 



K=0 



K=0 



We define the auxiliary operator p(s, t) as the Laplace trans- 
form of p(K, t) 



(8) 



K=0 



The long time behavior of the generating function Z has the 
following large deviation form 



Z(s,t) 



je(s) 



(9) 



The reason why this is the case will be more clear when we 
consider the dynamics of p(s,t) in the limit t — > 00. Equa- 
tion defines the function 8 which, in this formalism, is the 
analogue of a thermodynamic free-energy density, while s is 
interpreted as an intensive thermodynamic variable (the ana- 
logue of a chemical potential for example), conjugate to the 
counting variable K, and the time parameter plays the role of 
a volume. The function 6 is the fundamental object character- 
izing the dynamics of the process in terms of thermodynamic 
properties in the space of stochastic trajectories (the ensemble 
of possible time records for the counting variable K). Note 
that it is possible to consider not just one, but instead a set of 
jump processes associated to the jump operators {L ki }, and 
correspondingly a set of counting variables {Ki} can be taken 
into account. Accordingly one can introduce different inten- 
sive parameters {si}, conjugated to {Ki}. Since in the fol- 
lowing we associate one counting variable to each node of a 
graph, we use the vector notation for the variables s and K. 

The dynamical free-energy 8 can be obtained from the diag- 
onalization of a generalized Liouville superoperator describ- 
ing the dynamics of the auxiliary operator p(s, t). As shown 



3 



in appendix^ the Laplace transform in eq. dg)) decouples the 
evolution of p(K , t) for different vectors K, and this allows us 
to derive the following generalized quantum master equation 
for the auxiliary operator 

p(s, t) = Wfoft *)] = Hp& *)] + Ws, t)] . (10) 
L is the same Liouville superoperator as in eq. 



Nr 



h[.] = -i[H,.} 



k=l 



Lk-4-~{LlL k ,.}, (11) 



while 



v*H = 5>- d -i)£«i-4 
d 



d- 



(12) 



and d labels elements of a subset of {1, Nl}, those associ- 
ated to the jump processes that we want to track in time. When 
all the entries of s are we recover the real typical dynamics, 
otherwise W s - determines an evolution whose unfolding into 
trajectories is biased by e~ s . Following B32I1 we call this en- 
semble of trajectories the s-ensemble. In the classical context 
the analogue of Wj- is known as the Lebowitz-Spohn opera- 
tor lf34ll . and the free-energy function 8 is given by the largest 
real part of its eigenvalues. In the quantum case one can adopt 
a similar prescription and identify 8(s) with the largest real 
part of the eigenvalues of [32], accordingly with eq. @) 
and with the fact that for long times the largest real part of the 
eigenvalues will be the dominant one in evaluating Z(s,t). 
Once we have associated a dynamical free-energy to the en- 
semble of quantum trajectories of the system we can consider 
derivatives of the free-energy function, whose behavior tell 
us about the presence of dynamical crossovers or dynamical 
phase transitions. As an example calculation we now derive 
the form of the first and second derivatives of the dynamic 
free-energy. Using eqs. (0) to ®) we can write in the long 
time limit 

w = | LogIV(s - t) = MM. (1 „ 

The numerator in the last expression equals to 

Trp(s, t) = TrWzW, *)] = TrV s -[p(s, t)\, (14) 

where the last equality follows simply from eq. ^ and the 
cyclicity property of the trace. Hence we have that 

m = TlY 7 f^; t)] = Y,ie-"-V*Li d L d pfrt), (15) 
Z[s,t) - 

where we have introduced the normalized density matrix 

p(s,t) 



p(s, t) 



z(ity 



(16) 



Notice that 8(0) = 0, consistently with the existence of a 
steady state for the dynamics described by L. The partial 



derivative of 8 with respect to s d , for a particular d, is given 
by 

d d 8 = — ^ (e~ s * - l)Tr{Odd d p) - e - Sd Tr{O d p), (17) 

where O d = L d L d . We now introduce the activity vector 
parameter a, given by the inverse of the first derivative of 8(s) 



a d = -d d 8. 



(18) 



It is easy to see that, at s = 0, the activity coincides with 



Tr[O d p' 



Ml- 



/ ss ) being the stationary state associated with 



L. This means that, at s = 0, the dynamical activity a co- 
incides with the static expectation value of the operator O d in 
the steady state, so in this case there is a direct correspondence 
between static and dynamic order parameters. Hence it is not 
necessary to perform a dynamical analysis to access the value 
of the activity at s = 0. Nevertheless higher order derivatives 
of the free-energy function, and the behavior of the system 
at s ^ cannot be determined by static quantities. Those 
are the kind of information that can be accessed only with the 
LD approach. For a general value of s it is easy to show that 
ta d = [K d ] ave , in fact 



[Kd] 



tz 

-d d LogZ 



t 



d d t 



Otd- 



(19) 



This explains why a is named activity: it coincides with the 
typical number of jump events per unit of time, for a particular 
value of s. 

The second partial derivative of 8 in the d direction is given 

by 



d 2 d 6 - e- s «Tr(O d p) - 2e~ s «Tr(O d d d p) 
+(e- Sd - l)Tr(O d 8 2 d p). 



(20) 



For the second derivative, even at s = 0, one does not have 
a correspondence to any purely static observable. A useful 
quantity to look at in order to estimate the fluctuations in time 
of the local activity is provided by the index of dispersion 



6 d (s) 



_ dje(s) 

dd0(s) • 



(21) 



the ratio between the variance and the mean of the counting 
variable K for a particular jump process. In the next section 
we apply the above formalism in the context of dissipative 
continuous time quantum walks. 



IV. QUANTUM STOCHASTIC DYNAMICS ON DIRECTED 
GRAPHS 

In this section we study the dynamical properties of a QSW 
on a directed graph using the thermodynamic formalism. In 
particular, we consider the QSW that has been introduced in 
B35I1 with the purpose of defining a quantum analogue of the 
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classical notion of pagerank ll26ll . The latter is a particular 
graph centrality measure, i.e. an estimate of the relative im- 
portance of a node in a graph with respect to the dynamic of a 
random walker on the graph. Following [35] we consider the 
quantum master equation given by 

Nl , 

p = -i[H, p]+J2 l-.l'l-,, - r^'-.J-rl']- (22) 

fe=i 

where, in particular, H is the adjacency matrix of the graph 
obtained by removing the directions on the edges, and the 
jump operators are defined by 




U 



(23) 



The matrix G, called Google matrix, is a stochastic matrix as- 
sociated to an irreducible and aperiodic classical random walk 
on the graph (for a precise definition of the Google matrix see 
l27lo . As shown in 13511 . using a theorem by Spohn II36J it 
is possible to prove that the dynamics described by eq. (1221) 
is relaxing, i.e. that the stationary state is unique. The clas- 
sical random walk described by the stochastic matrix G also 
has the same property, and its stationary state is called pager- 
ank, denoted with 7?. Notice that, without the coherent part 
given by H, the dynamics for the diagonal elements of p(t) 
provided by eq. (1221) is identical to the classical random walk 
described by the Google matrix, and in particular its stationary 
state coincides with the classical pagerank 13511 . On the other 
hand, adding the coherent term H modifies the dynamics and 
the stationary state, allowing for a quantum generalization of 
the pagerank vector. We are now going to use the thermody- 
namic formalism to detail the differences between the purely 
classical and the more general quantum open dynamics. 

We start calculating the value of the activity at s = for 
the set of jump operators {Ljj; j — 1, n}, associated to a 
chosen node i. Using eqs. (U7I) and (1231) we have 



a i (0) = Y / G ij Tr[\j){j\p^) = ^2G 



(ss) 

IJPjj ■ 



(24) 



This equation tells us that the activity at s = is related 
in a simple way to the diagonal entries of the steady state, 
through the stochastic matrix G. Without the presence of the 
Hamiltonian H in eq. d22h we would have that p( ss ^ is equal 
to diag(7Ti, . . . , 7r„), hence 



<*(0) »</C 



(25) 



where the last equality follows from the definition of the 
pagerank vector, the stationary state of G. This result is sim- 
ply a check of the consistency of the thermodynamic formal- 
ism applied to the classical pagerank, and it states that at long 
times the probability of finding the walker in the z-th node is 
equal to the average number of times the walker hits the i-th 
node. It is interesting to observe that the quantum case is dif- 
ferent. When the coherent term H is present the stationary 
state is different from the mixed classical state associated to 



FIG. 1. The directed graph that we consider as an example for the 
application of the thermodynamic formalism to QSWs. 



the pagerank vector: p^^ ^ diag(7Ti, . . . , 7T n ). Hence the 
equation 



(ss) 
33 



(26) 



defines a dynamical quantum analogue of the classical pager- 
ank given by a, which is different from the static quantum 
pagerank p( ss ) of lf35ll . The values provided by a quantify 
the typical number of times the walker jumps into a particu- 
lar node in the long time limit. On the other hand p( ss ^ pro- 
vides the probability of finding the walker at a particular node. 
While a purely classical dynamics would not distinguish be- 
tween the two quantities, as shown by eq. J20I) . the quantum 
dissipative dynamics does differentiate the two because of the 
coherences (off-diagonal elements) present in the steady state. 

We consider now, as an example, the QSW on the directed 
graph shown in fig. [T]. We set the dimension of the vector s to 
be equal to the size of the graph and, in order to simplify the 
discussion, we choose all the entries of s equal to the same 
parameter s. A similar analysis for a simpler 2-node directed 
graph can be found in appendix[l|. The dynamical free-energy 
9 is plotted in fig. where one can see that it is equal to 
at s = 0, consistently with the fact that the physical dynamic 
is rela xing , given that 8 is the largest real eigenvalue of W s 
in eq. (flOb . More interesting is the behavior of the first two 
partial derivatives of the free-energy function. In fig. H] we plot 
the normalized activity with respect to different sites 



(n) _ 0!i(s) 
H = 6 ' 

Ei=i a i(«) 



(27) 



The figure shows a smooth dynamical crossover from the ac- 
tive region, where large negative values of s favour trajectories 
with many jumps, to the inactive region, where large positive 
values of s favour trajectories with few jumps [see eq. (fl2l)l. 
Note that eq. d27l) is not well defined in the limit s —> 00, since 
in this case one is left with a purely deterministic evolution ac- 
cording to H e g [see eq. (Q)], which implies cti(s — > 00) = 0. 



-3-2-1 1 2 3 

S 

FIG. 2. The dynamical free-energy 8(s) for the QSW associated to 
the graph in fig.Ql 



FIG. 5. (color online) The stability parameter quantifying fluctua- 
tions in the trajectories for different values of s. Each line corre- 
sponds to a node in the graph of fig.Q]. The inset shows the sum with 
respect to all nodes of the same quantities, and it shows a maximum 
close to s — — 1. 




FIG. 3. (color online) The normalized dynamical activity of the six 
nodes of the graph in fig. Q]. For this particular example the rank- 
ing provided by the activity order parameter does not change as a 
function of s (see text for more details). 
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FIG. 4. (color online) A comparison between the classical pagerank 
(left bars), the dynamical activity at s — (middle bars), and the 
populations in the steady state (right bars). These values are obtained 
with respect to the graph shown in fig. Q]. Note that the extreme left 
values in the plot correspond to the classical pagerank. 



This is not a problem for our approach, since one could al- 
ways use the unnormalized activity a, which is well defined 
for every value of s. However the normalized is useful, 
since it makes the quantitative comparison between the classi- 
cal pagerank and its quantum dynamical generalization easier. 
Moreover, in the limit s — > — oc the dynamics in eq. (TTob is 
dominated by with no time intervals of deterministic evo- 
lution. These two extreme scenarios select trajectories that 
are either associated to "lazy" QS Ws which never jumps, or 
to very active QSWs that are effectively equivalent to the clas- 
sical random walk described by G (in fact in this limit co- 
incides with the classical pagerank). The transition between 
the two extreme ensembles of trajectories is controlled by the 
parameter s. For the particular example that we consider the 
activity ranking (the sorted list of the elements in the vector 
a) is preserved for all values of s, so in particular it coincides 
with the classical ranking provided by n, as shown in fig. |H 
However the specific values of the activity ranking is a func- 
tion of s, and it could be the case that particular graphs show 
some rank crossings as s is changed. We note that the situ- 
ation is different if one considers the quantum ranking pro- 
vided by the population (the occupation probability) at each 
node in the steady state lf35ll . As shown in fig. the popula- 
tion ranking is different for nodes 2 and 3, with respect to the 
classical pagerank and the activity ranking. This simple exam- 
ple shows a peculiar aspect of QSWs with respect to classical 
random walks, i.e. their richer dynamical behavior translates 
in distinguishable ways of estimating nodes' centralities. 

We now look at the dynamical fluctuations of the counting 
processes. In fig.dwe plot Si [see eq. J2ll) l with respect to dif- 
ferent nodes. As we said before Si is a measure of the disper- 
sion of the counting variable with respect to its average value, 
and it provides information on the fluctuations in time of the 
stochastic process. In particular for Si > 1 the process is said 
to be over-dispersed, meaning that the jumps are clustered in 
time, having time windows with many jumps, and time win- 
dows with almost no jumps. On the other hand, if Si < 1 the 
process is under-dispersed, meaning that the jump events are 
more uniformly distributed in time. The dynamical index of 
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FIG. 6. (color online) The figure shows the activity order param- 
eter for the 4-th node in the graph of fig. Q] (triangles), and the total 
fluctuations of the dynamical trajectories on the same graph (circles). 
The presence of a maximum in the fluctuations corresponding to the 
region where the dynamical order parameter is changing is an indi- 
cation of a crossover phenomenon. 



dispersion Si can itself provide a centrality measure for the i- 
th node, quantifying the stability (or robustness) of the activity 
index with respect to fluctuations: the smaller Si, the smaller 
are the fluctuations in time of a,i. Similarly one can introduce 
a global index of dispersion: 5 = Si, quantifying the fluc- 
tuations at the level of the whole network. In fig. y one can 
see a crossover in the dynamical robustness from the active 
region (associated to an effectively classical dynamics), to the 
inactive region. Comparing the two extreme behaviours one 
notices a remarkable difference in terms of dynamical fluctu- 
ations of the local activities. In the classical regime the top 
3 ranked nodes are also the more stable, while the lowest 3 
are over-dispersed. In the quantum regime this is not the case 
anymore, and all the nodes fluctuates with more similar in- 
tensities. In general their stability ranking is strongly depen- 
dent on s in a region approximately delimited by s = — 2 and 
s = 0. The inset in the same figure shows the total index of 
dispersion, i.e. the sum of all site's indexes, which measures 
the global dynamical fluctuations. It is interesting to observe 
that in the transition from classical trajectories to quantum tra- 
jectories the behavior of the total dispersion is not monotonic, 
and it shows a maximum close to s = — 1 in the active region. 
The maximum in the global index of dispersion supports the 
interpretation of a thermodynamic crossover from an effec- 
tively classical to a quantum regime in the space of trajecto- 
ries. This can be seen more clearly in fig. |(| where we plot the 
relation between the drop in the normalized activity of the 4th 



sites (the most active), and the global dispersion index. From 
the figure one can identify two distinct values of the activ- 
ity order parameter, corresponding to distinguished dynamical 
phases, together with an increase in dynamical fluctuations in 
the region separating the two, which includes the physically 
relevant point s = 0. In fig. d one can see that the stabil- 
ity ranking is highly dependent on s around 0, which implies 
that small fluctuations in the dynamical trajectories will have 
a noticeable effect on Si, 

V. CONCLUSIONS AND OUTLOOK 

In conclusion, we have applied the large deviation approach 
to study the dynamics of quantum stochastic walks, a particu- 
lar case of dissipative quantum evolution on directed graphs. 
The LD approach allows one to formulate a thermodynamic 
formalism describing quantum dissipative evolutions. We find 
that the properties of the ensemble of dynamical trajectories 
of the quantum walker can be used to characterize the struc- 
ture of the underlying graph. We introduce two new quan- 
tum centrality measures 0711 . the activity rank and the stabil- 
ity rank. These measures are of relevance in the context of 
complex networks theory ll38l - l40ll . providing another connec- 
tion between complex networks and quantum information the- 
ory ll35i I4ll - l55ll . Another interesting aspect of the thermody- 
namic approach to dissipative quantum walks is the possibil- 
ity to frame dynamical phases of the evolution in the standard 
language of Statistical Mechanics and Thermodynamics. We 
provide an example showing that, already on a small directed 
graph, the quantum walker displays two dynamical regimes 
separated by a smooth thermodynamic crossover. Note that 
this approach is valid also for undirected graphs, though in 
this case the identification of the most relevant nodes can be 
simply done looking at the degree distribution: the higher the 
degree of a node, the higher will be its classical or quantum 
centrality measure; while the information related to fluctua- 
tions around the typical values is less trivial, and it can be 
obtained with the LD approach. 

It would be of interest to understand the connections be- 
tween the topology of the underlying graph and the dynamical 
phase space of the quantum dissipative process. The compar- 
ison between the simple case of a two-node graph provided 
in appendix |b| -where no crossover is present- and the more 
structured graph in fig. [II suggests that the topology of the 
underlying graph can be responsible for the occurrence of dif- 
ferent dynamical phases. In particular, exploring more com- 
plex graph structures, one could characterize the cases where 
the transition is sharp (as found in other physical contexts 
dU-Ult]), which would imply a dynamical phase transition be- 
tween effectively classical and quantum coherent regimes. 



[1] T. L. Hill and Physics, An Introduction to Statistical Thermody- 
namics (Dover Books on Physics) (Dover Publications, 1987). 

[2] S. Sachdev, Quantum Phase Transitions, 2nd ed. (Cambridge 
University Press, 201 1). 

[3] F. Verstraete, M. Wolf, and J. Cirac, Nature Physics 5, 633 



(2009). 

[4] S. Diehl, A. Micheli, A. Kantian, B. Kraus, H. Buchler, and 

P. Zoller, Nature Physics 4, 878 (2008). 
[5] B. Kraus, H. Buchler, S. Diehl, A. Kantian, A. Micheli, and 

P. Zoller, Physical Review A 78, 042307 (2008). 



7 



[6] S.Diehl, A. Tomadin, A. Micheli, R.Fazio, and P. Zoller, Phys- 
ical review letters 105, 15702 (2010). 

[7] D. Aharonov, A. Ambainis, J. Kempe, and U. Vazirani, in Pro- 
ceedings of the thirty-third annual ACM symposium on Theory 
of computing (ACM, 2001) pp. 50-59. 

[8] E. Farhi and S. Gutmann, Phys. Rev. A 58, 915 (1998). 

[9] A. Childs, E. Farhi, and S. Gutmann, Quantum Information 
Processing 1, 35 (2002). 
[10] E. Farhi, J. Goldstone, and S. Gutmann, Theory OF Computing 
4, 169 (2008). 

[1 1] A. M. Childs, Phys. Rev. Lett. 102, 180501 (2009). 

[12] O. Mulken and A. Blumen, Physics Reports (2011). 

[13] D. I. Tsomokos, Phys. Rev. A 83, 052315 (2011). 

[14] P. Rebentrost, M. Mohseni, and A. Aspuru-Guzik, The Journal 

of Physical Chemistry B 113, 9942 (2009). 
[15] M. Mohseni, P. Rebentrost, S. Lloyd, and A. Aspuru-Guzik, 

The Journal of Chemical Physics 129, 174106 (2008). 
[16] M. Plenio and S. Huelga, New Journal of Physics 10, 113019 

(2008). 

[17] F. Caruso, A. Chin, A. Datta, S. Huelga, and M. Plenio, The 
Journal of Chemical Physics 131, 105106 (2009). 

[18] P. Rebentrost, M. Mohseni, I. Kassal, S. Lloyd, and A. Aspuru- 
Guzik, New Journal of Physics 11, 033003 (2009). 

[19] P. Giorda, S. Garnerone, P. Zanardi, and S. Lloyd, Arxiv 
preprint arXiv: 1106. 1986 (2011). 

[20] F. W. Strauch, Phys. Rev. A 79, 032319 (2009). 

[21] F. Strauch and C. Williams, Physical Review B 78, 094516 
(2008). 

[22] J. Whitfield, C. Rodriguez-Rosario, and A. Aspuru-Guzik, 
Physical Review A 81, 022323 (2010). 

[23] M. Plenio and P. Knight, Reviews of Modern Physics 70 (1997). 

[24] H. Touchette, Physics Reports 478, 1 (2009). 

[25] D. Ruelle, Thermodynamic formalism: the mathematical struc- 
tures of equilibrium statistical mechanics (Cambridge Univ Pr, 
2004). 

[26] S. Brin and L. Page, Computer Networks and ISDN Systems 
30, 107 (1998). 

[27] A. N. Langville and C. D. Meyer, Google's PageRank and Be- 
yond: The Science of Search Engine Rankings (Princeton Uni- 
versity Press, 2006). 

[28] Y. Shikano, in AIP Conference Proceedings, Vol. 1327 (2011) 
p. 487. 

[29] S. Venegas-Andraca, Synthesis Lectures on Quantum Comput- 
ing 1, 1 (2008). 

[30] H. Breuer and F. Petruccione, The theory of open quantum sys- 
tems (Oxford: Oxford University Press, 2002). 

[31] L. Hedges, R. Jack, J. Garrahan, and D. Chandler, Science 323, 
1309 (2009). 

[32] J. Garrahan and I. Lesanovsky, Physical review letters 104, 
160601 (2010). 

[33] J. Garrahan, R. Jack, V. Lecomte, E. Pitard, K. van Duijvendijk, 
and F. van Wijland, Physical review letters 98, 195702 (2007). 

[34] J. Lebowitz and H. Spohn, Journal of Statistical Physics 95, 333 
(1999). 

[35] E. Sanchez-Burillo, J. Duch, J. Gomez-Gardenes, and 

D. Zueco, Scientific Reports 2, 605 (2012). 
[36] H. Spohn, Letters in Mathematical Physics 2, 33 (1977). 
[37] S. Boccaletti, V. Latora, Y. Moreno, M. Chavez, and D.Hwang, 

Physics reports 424, 175 (2006). 
[38] A. Barrat, M. Barthlemy, and A. Vespignani, Dynamical 

processes on complex networks (Cambridge University Press, 

2008). 

[39] A.-L. Barabasi, Linked: The New Science of Networks, 1st ed. 
(Basic Books, 2002). 



[40] A. Barabasi, Nature Physics 1, 68 (2005). 

[41] O. Giraud, B. Georgeot, and D. Shepelyansky, Physical Review 
E 72, 036203 (2005). 

[42] S. Braunstein, S. Ghosh, and S. Severini, Annals of Combina- 
torics 10, 291 (2006). 

[43] H. Kimble, Nature 453, 1023 (2008). 

[44] G. J. Lapeyre, J. Wehr, and M. Lewenstein, Phys. Rev. A 79, 
042324 (2009). 

[45] S. Perseguers, M. Lewenstein, A. Acm, and J. Cirac, Nature 

Physics 6, 539 (2010). 
[46] S. Perseguers, Phys. Rev. A 81, 012310 (2010). 
[47] S. Perseguers, D. Cavalcanti, G. J. Lapeyre, M. Lewenstein, 

and A. Acm, Phys. Rev. A 81, 032327 (2010). 
[48] K. Anand, G. Bianconi, and S. Severini, Phys. Rev. E 83, 

036109 (2011). 
[49] L. Wu and S. Zhu, Phys. Rev. A 84, 052304 (201 1). 
[50] M. Cuquet and J. Calsamiglia, Phys. Rev. A 83, 032319 (201 1). 
[51] R. Ionicioiu and T. P. Spiller, Phys. Rev. A 85, 062313 (2012)] 
[52] S. Garnerone, P. Zanardi, and D. A. Lidar, 

|Phys. Rev. Lett. 108, 230506 (2012)) 
[53] G. D. Paparo and M. A. Martin-Delgado, Scientific Reports 2, 

444 (2012). 

[54] S. Garnerone, P. Giorda, and P. Zanardi, New Journal of 

Physics 14, 013011 (2012). 
[55] B. Adhikari, S. Adhikari, and S. Banerjee, Arxiv preprint 

arXiv: 1205.2747 (2012). 



Appendix A: Derivation of the generalized quantum master 
equation 

We denote with 

K=[Ki,...,Ki,...,Ki\ (Al) 

the vector corresponding to the occurrence, up to a certain 
time, of a number of Ki events associated to the jump pro- 
cesses described by Lf. K\ jumps associated to L\, K% jumps 
associated to L-2 , and so on. The probability at time t for the 
set of events K to occur is given by 

p(K;t) = Trp(K;t), (A2) 

and p(K\ t) is the reduced density matrix of the system con- 
ditioned on the observation of the set K of events. We define 
the generating function Z as follows 

oo 

Z(s;t)= e-* it p(K;t)=T I p{s;t), (A3) 

K u ...,K l= 

where 

p(?;f) = £y**p(J?;f). (A4) 

it 

Using eq. (Q]) the time evolution of p(K;t), in the K sub- 
manifold, is given by 

p(K; t) = L [p(K: <)] + ]T L, ; *)] , (A5) 
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with 

ho['] = -i[H,-]-^{LtL ir }, 

i 

U[p{K (l) ;t)]=Y,L l p(K 1 ,...,K l -l,...,K l ; t)L\ . 

i 

(A6) 

It follows that the generalized quantum master equation for 
p(s; t) is given by 

p{s;t) = Y t e- 3 - R p{K ] t) = 
it 

E e~**U [p(K; t)} + E e~** E L < \P$® > *)] = 



K 



L 



Y,e- SK p{K-t) 
it 



E L > 



5>-**p(#0;t) 



K 



e"^^' ]T e- s ^-Vp(K^;t) 



xL 



Ki=X 



E e-^^Y.z- S ^p{K,t) 

{Kj^i}=0 K i= Q 

L [p(s; t)] + E e~'*U [p(s; t)} = W a - [p(s; t)] . (A7) 



In the above derivation we simply made use of the linearity of 
the operators, and in the third equality we wrote the YIk =0 
into a second equivalent way. 



shows the result of the simulation. The function 9 has the ex- 
pected monotonic decreasing behavior as a function of s. The 

(n) 

normalized activities a 2 of the two nodes deviate most go- 
ing from the active to the inactive dynamical region, which is 
different from the case of fig. [H where more active trajec- 
tories allow to better distinguish the activities of the nodes. 
The index of dispersion <5,-i.2 increases in a smooth way for 
both nodes from the active to the inactive region, approach- 
ing a limiting value of 1, characteristic of a Poisson process. 
The global index of dispersion S does not show any maximum 
between the active and the inactive region, implying the ab- 
sence of dynamical crossovers between the two regimes, sup- 
porting the conjecture that more interesting dynamical effects 



O — 



FIG. 7. A directed two node graph. 




Appendix B: Two-node directed graph 

The most simple example of a directed graph is given by a 
two-node graph with a directed link (see fig. 0). This example 
is useful because it shows that a trivial graph structure does 
not show the dynamical crossover present in fig. [(§ Figure |g] 



FIG. 8. (color online) Simulation results for the directed 2-node 
graph in fig.0. See the main text for comments. 



(like crossovers or phase transitions) are associated to more 
structured graphs. 



